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Abstract 

We show that Liouville-von Neumann approach to quantum mechanical sys¬ 
tems, which demands the existence of invariant operators, reproduces the time- 
dependent variational Gaussian approximation. We find the effective action of the 
time-dependent systems and show that many aspects of the dynamics are inde¬ 
pendent of the details of time evolution, e.g., the squeezing of the wave-function 
is determined by the effective potential of the final stage of time-evolution. 
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1 Introduction 


The interests in searching for a method which effectively describes the real-time quantum 
dynamics has been continually increased, in relation with the quantum theory out of 
thermal equilibrium [Q], the inflationary scenario of the Universe 0, the dynamic process 
in phase transitions |P], the formation of topological defects in variety of phenomena in 
condensed matter physics Q]. Understanding quantum mechanical evolution of inflaton 
held in the new inhationary universe is vital for correct determination of the inhationary 
epoch’s period [Q. 

A time-dependent variational approximation for quantum dynamics was developed 
in Schrodinger picture |^. The Liouville-von Neumann (LvN) approach of the time- 
dependent harmonic 0] and driven harmonic 0 oscillators have been developed to obtain 
the exact time-dependence of quantum states, and the methods are generalized to the 
cases of anharmonic oscillator, scalar inhaton held dynamics |10], and driven anharmonic 


oscillator 0. LvN approach has a merit that it can be directly applied to both time- 
independent and time-dependent quantum systems without any modihcation, because 
it is based on the operator level equations. 

The purpose of this paper is to extend the approach developed in Ref. to include 
arbitrary potential, so that one can develop a general Gaussian approximation using the 
LvN method. This approach has two merits: 1) it can describe the time-evolution of 
systems from the start to the end. 2) In LvN approach, it is easy to dehne the general 
Fock space because it directly identihes the creation and annihilation operators. 

The organization of this paper is as follows. In Sec. II, the LvN approach is elabo¬ 
rated to be applicable to the case of general potential V(q, t) in obtaining the Gaussian 
approximation, and is applied to hnd the Fock space and the time-dependent expecta¬ 
tion values of some operators for V{q,t). Sec. Ill is devoted to the study of a model 
with slow rolling potential. In Sec. IV, we summarize our results and present some 
discussions. 


2 Gaussian approximation in Liouville-von Neumann 
approach 

The quantization of a time-dependent anharmonic oscillator using the LvN approach 
was studied in Ref. [^. It was shown that their approach for symmetric potential is 
equivalent to the variational Gaussian approximation and the mean-held, Hartree-Fock 
method. In this section, we develop LvN method for a general potential and show that 
this approach exactly reproduces the Gaussian approximation of Ref. [0, without losing 
the merits of the LvN approach. 

We start with the Hamiltonian: 


Hit) = 


2m(t) 


+ V{q{t),t), 


( 1 ) 
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where the potential is an arbitrary function of q{t) and t. 

Now let us assume the time-evolution of an initial ground state can be described by 

2 ' 

a Gaussian wave-packet which is annihilated by an operator: 

A = i[u*{t)p{t) + v*{t)q{t) + X*] . (2) 


The LvN approach is to demand this operator to be invariant under time evolution: 


dA 

dt 


• * I 

u + 


V* 

m{t) 


p{t) + i [v*q{t) 


u*V\q,t) + X*] . 


( 3 ) 


Taking the expectation values with respect to the wave-function T on each side of the 
equations obtained from the variations in p or g of Eq. (^, gives 


v*{t) =-m{t)u{t), v{t) -u{t){V''{q,t)) = 0, ( 4 ) 

x*{t)+u*{t) [{q{t)){V”{q,t)) - {V'{q,t))] = 0 . ( 5 ) 


These equations reduce to one equation for u: 

m{t)u+ {V" {q,t))u{t) = 0 ) 


( 6 ) 


where u(t) = ^Jm(t)u(t). This equation is too complexly coupled to x(t) to be solved 
exactly, except for a simple problem of harmonic oscillators 0. To overcome this 
difficulty and to hnd a way for easier application to physical problems, we dehne an 
auxiliary function 


q{t) = —ih [u{t)x*{t) — u*{t)x{t )], 
which satishes the equation of motion 

m{t) 


( 7 ) 


( 8 ) 


The equation (|^) becomes the ordinary differential equation for q only in the semi- 
classical approximation. In general, the expectation value {F{q,t)) depends on q and 
= hlhp. Therefore, we need also to consider the equation of motion for r 




r 


m{t) 


( 9 ) 


where L = r‘^{t = to)4>{t = to) is a constant of motion, with 0 being the phase part 
of u{t). If the initial wave-packet at origin is a static Gaussian, the initial condition 
becomes Ui(t) = and a:(0) = 0, which gives L = —1/2. The two equations 

(H), (H) and the boundary condition mentioned above determine a unique solution for q 
and r. 
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Once we get u and x from Eqs. (|^), (^, (H), and (|^), the time dependence of quantum 
dynamical operators is explicitly written as: 


m{t) 


q(t) = u{t)A-\-u*(t)A^+ 


p(t) = m{t) \u{t)A +u*(t)A^ + 


( 10 ) 


d q{t) 


Therefore, the expectation values of operators at time t for the state T are given by 

(11) 


(■ill)) = 


m{t) 


dt 


m{t) 


m)) = 

Qit) 

2 

+ hM^ {p^(t)) = 

_Jm(t)_ 




m{t) 


d q{t) 


dt 


1 2 


m{t) 


+ {t)h^\u\‘^. ( 12 ) 


To have explicit form of Eqs. (|^) and (H), we need to expand the expectation value of 
operator F{q, t) in series around q: 






n=0 


n\ 


q{t) - q/\]m{t) 


where the expectation values of the powers of q{t) — q/qjm{t) 
the iteration formula: 


). (13) 

can be obtained from 


( 


q{t) - q/^Jm{t) 




q{t) -q/^Jm{t) 


n—2 


), n>2. (14) 


Note that this series is the same as the series expansion in h because h'^\u\‘^ is 0(h). 
The equations (^ and (^ become to 0(h), 


x*(t) + u*(t) 


^ ^(2) _ yiO y f V"(h _ yiA 


m 


2m \ \/m 


= 0 , 


(15) 


m(t)u(t) + 


2 l -^|2 


y(2) ^ ^ \U\ .^^( 4 ) 

2m 


u(t) = 0, (16) 


where = d!^V(x,t)\,^ =qlA>n independent of r. The induced set of equations 

and (H) is given by 


-i,(t)+A==+o^v'^^ = 0 . 




(17) 

( 18 ) 
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These equations of motion admit to define the effective Hamiltonian, 


Heff{q,r) = 


^2 

q 


2 2 2 8f2 8 

j + y +W9,i)) + 5^. 


( 19 ) 


where f = = h|M|. Eq. ([T^), which is a classical Hamiltonian for an anharmonic 

oscillator in two dimensional space x® = (g, f), represents the same effective potential 
of Cooper, Pi, and Stancioff 0]. In this sense, the present approach exactly reproduces 
the variational gaussian approximation of previous authors. The LvN method, however, 
presents better understanding of the quantum mechanical system, because not only it 
naturally dehnes the complete Fock space through A and but also these operators 
are time independent. In this sense, the parameterizations of the wave-function by u 
and X is more convenient than by q and r. 

From the vacuum dehned by HlO) = 0, we can naturally dehne the excited states by 
successively applying the creation operator on the ground state. The mixed states 
can also be dehned. For example, the thermal states can be dehned from the density 
matrix [ITU 


p^(t) = — 

Zn 


( 20 ) 


where N{t) = (3 is the inverse temperature, ujq is the frequency of u at f = 0, and 

the partition function Ztv is given by 


OO 

n=0 


1 

2sinh(^)' 


( 21 ) 


We explicitly write down the expectation values of some operators with respect to the 
thermal state: 


{q{t))T = g(t), (5g^(t))r = ^V(^)Pcoth^^, (22) 

{p{t))T = (^p^(t))r = h^?^^|h(t)|^ coth 

3 Slow rollover transition 


We set m{t) = 1 in this section, since the mass dependence can be removed by reparametriza- 
tion of other parameters in the Hamiltonian. We thus consider the potential 

V{qA) = {[<?W - At)f + Ht)} , (23) 


4 










where we assume a^(t) [> —A;(t)] increases from zero to some constant number a^, and 
k{t) decreases from some positive number to a fixed non-positive number k at t — >■ cx), 
which assures the potential at g = 0 to be locally stable false vacuum. The dynamics for 
this potential contains many interesting quantum mechanical features depending on the 
time-dependence of k{t) and a{t), such as tunnelling (—< k < 0), phase transition 
(k = 0), or slow rollover (k = —a^). 

The energy for a static Gaussian wave-packet at f = 0 described earlier is given by 


E — —hujQ -\- -\- ki) -\- 

4 48 


32a;^‘ 


(24) 


The effective potential corresponding to Eq. 


yeff{q,r) = V{q,t)+ '-^ q^ - aq + 


T 


becomes 
+ k 


6 


^ 8f2 


(25) 


In the case of a time-independent system, we can describe the whole dynamics in terms 
of this effective potential. We present a contour plot of the effective potential in Fig. 1. 



Fig. 1. The equipotential curves of the effective potential. The horizontal axis represents q and 
the vertical axis r. In this figure we use the parameters, a = 8, k = —10, h^/X = 0.1, and each curve 
represents the energy values, from the innermost one to outer direction, U^E/X =0, 1, 2, 3, 4, 5. A 
given Gaussian wave-function cannot escape the contour defined by its energy. 
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One crucial difference of the present paper from that of Cooper et. ah is that we 
keep f dependence on the effective potential, while it was variationally determined in 
Ref. 0. This results in a large difference in the slow rolling period. As clearly shown 
in Fig. 1, a given wave-packet has allowed region of {q,f) determined by the value of 
its energy. Depending on the initial value of the wave-function simply oscillates 

around the false vacuum g = 0 or go over to the true vacuum near q = 8. There are 
two properties which are distinguished from its classical dynamics. First, there is a 
quantum mechanical driving force at g = 0. The classical stable equilibrium at g = 0 
is not the equilibrium any more in quantum mechanics. Second, it was shown that 
as the dispersion (frequency) becomes smaller (larger), it becomes easier for g to go 
over the classical potential barrier |p. In Fig. 1, it is evident that as the dispersion 
(frequency) becomes larger (smaller), it also becomes easier for g to go over the classical 
potential barrier. There is an inevitable restriction of the Gaussian approximation. In 
the presence of leaking into the right well the wave-function cannot be described by the 
Gaussian. Due to this restriction, we cannot effectively describe tunnelling process. To 
remove this restriction, one needs double Gaussian (Appendix B of Ref. [^) or other 
approximations, which is out of scope of the present paper. If a{t) becomes time- 
dependent, the dynamics cannot be understood simply by the effective potential, and 
the equations of motion should be considered explicitly. 

As an explicit dynamical example, we consider the slow rolling {k + a"^ = 0 for t > 0) 
potential. The equipotential curves are shown in Fig. 2. An interesting features of this 
effective potential is that the unstable equilibrium point at g = 0 disappears. 

To understand the dynamics of this system, we write down the equation of motion 
for X* = (g, f) explicitly: 


X 


+ = 0 , 


(26) 


where Veff{x^) is given by Eq. (p5|). Let a Gaussian wave-function prepared at g = 0 
with velocity g(0) = 0 starts to evolve. In the case of slow rollover, the evolution is 
divided into three phases; the slow-rolling, rolling, and oscillating phases. 

As soon as the packet released, it starts to roll down the potential wall. In the region 
g ~ 0, f increases slowly until f ~ [h^/ (2A)]^/®, which is the variational minimum of the 
potential in f at g = 0, because f ~ ^ (2ix;o)^^^ and f(0) = 0. This time (ti), which 

f{t) takes in reaching to its variational value, is roughly given by ti ~ ^ ' 

After that, the time evolution of (g, f) is similar to its variational path. The initial time 
evolution of g is governed by the quantum-mechanical driving force given by a{t) until 
(f ~ f^. Therefore the time (^ 2 ) for g to increase to f, can be obtained from 


g(G) - \j ~ r{t2) ~ r(g). 


(27) 
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Fig. 2. The equipotential lines of a slow rolling transition. The horizontal axis represents q and the 
vertical axis f. In this figure, o = 6, ?i^/A = 0.1, and each curve represents the energy values 
Ti^E/X =10, 0, -10, -20, -30, -40 from outside. A given Gaussian wave-function cannot escape the 
contour defined by its energy. 

If ti > t 2 , we determine t 2 from Eq. (|^) , which gives t 2 ~ ^ • Otherwise, 

t 2 is given by the sum of the time taken by f{t) to increase to the variational value, ti, 
and the variational time for the remaining part of q, which will be shown shortly. 
The variational effective potential is Kar = K//(?,h(g)), where f(g) satishes 

f®(g) -1- if - af)f^{q) - ^ = 0. (28) 

Therefore, in the variational approach, g(0) = ^ , and the duration of slow- 

rolling phase is ~ Note that the duration (ts = ti or tg = t 2 +1^) of the slow- 

rolling phase in both of the scenario is considerably larger than the time determined by 
the variational method 

For t > tg, the wave-packet runs into the rolling phase. In this phase, the equation 
of motion becomes too difficult to get exact answer without numerical analysis. In this 
paper we simply outline the evolution without giving explicit solution. If a{t) = a is a 
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constant for t > tg, and is smaller than q^, one gets q{t) from the elliptic integral 


dq A , , 

, = —(t - tg), 

sjc + q^{-2aq + q^) 12 


( 29 ) 


• 2 

where c = 12q {tg)/\. However, as one can see from Eq. (|^) , f grows qnickly dnring 
the period < aq — q^, which lasts nntil q ' 2 :^ a. Once f or g grows enongh and fails to 
satishes the condition, < aq — q"^, f goes into the oscillating phase. 

After q reaches to the potential minima, it starts to oscillate aronnd the minimnm of 
the potential, which is the last stage. We do not restrict to the slow rolling case in the 
calcnlation of this stage. We simply assnme that the hnal wave-packet resides aronnd the 
trne vacunm, which is determined by VW//|so = 0. For a slow rolling potential, there 
is a nnique point which satishes this condition. Now, we series expand the potential 
aronnd the trne vacunm: 

Vt{x) = Veffixo) + ^ 51 didjVeffix)\^=a,o{x' - xl){x^ - x^) + ■ ■ ■. (30) 

The solutions for the equation of motion for this Vr are given by two frequency series 
expansion in general: 


X\t) = Xl[l + '^ Xn,m cos n{uJmt + 0m)], 

n=l m=± 


(31) 


where Xn,m is determined from the differential equation (PB]) up to its scale, and uj± is 
the normal mode frequencies of Vt- We show this fact to the hrst order in the series 
expansion. Using the coordinates transformation: 


= 


q + br 


X_ = 


bq — r 


VTT¥' ViT¥' 

the Hamiltonian, up to quadratic order, becomes 

Ht = xI + xI + - Lixi + 


where b and uj± satisfy 
b^ - 


2 L 


2 dq^eff{x)\x=xo dr^ef f{x)\x=XQ 


2did2Veff{x)\x=, 


6-1 = 0 , 


XQ 


= 




eff\^)\x=XQ 


1 + 62 


^2 (^) I x=xo ■ 


(32) 


(33) 


(34) 


(36) 


The solutions for the Hamiltonian system (|33[) are X±(t) = Xo,± cos(a;±t + 6*±), where 
Xo,± should be determined from the whole time evolution. From this we get, to hrst 
order in Xq, 

XJt) + bX_(t) bXXt)-X(t) , , 

q(t)=qo + ^^==^, f(t)=ro+ (36) 


VTT¥ 


VTT¥ 
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We have two unknown parameters Xq^± which should be determined from the whole time- 
evolution. If the system is time independent for f > 0, one of them can be determined 
from the energy conservation law. 

We now turn our attention to the calculation of u{t) = This can be done 

by integrating 0 = l/(2r^). But, this integration is highly nontrivial if Xq is not much 
smaller than ro, and the slow rolling is the case. As a round about method, we directly 
try to solve the equation for u, and then determine the free parameters in u from the 
effective potential. To simplify the procedures, we assume 6 -C 1 and ro -C q'o- These 
two assumptions are usually satished for large %. Now we write 


« = 9o[l - xmsiujt + 9)], 


(37) 


where qox 
to be: 


—Xo,+. The equation ([T6|) for u, to linear order in y, can be approximated 


u{t) -f- a;^[l — x' cos(a;t -|- 0)]u{t) = 0, 


(38) 


A(9o - C‘qo/2) 




X- With the transformation ujt + 9 = 20, Eq. (^) becomes 


where x' = 

the Mathieu equation with s = 2, 6 = 4 -|- 4y', and h? = 8y' in page 562 of Ref. [|T2 
Its solution, to the hrst order in y', is given by 


u{t) = auj^it) + 


(39) 


where a and B are constants, and u+it) = , 

that M+ satishes Uj^u\ — u\uj^ = i. Note also that 

f^(t) = h^\u{t)\^ = h^[(|Q(|^ -1- \(5\^)\u+{t)\^ a(3u\ -\- a* (40) 

where its time average h{\a\^ + \j3\‘^)/{2uj) should give fg. From this relation and Eq. (^21) 
below, we determine \(3\ = y^curg/h — 1/2. Yet undetermined parameters are y, the size 
of fluctuation of q, and the relative phase between a and 0, whose determinations need 
the analysis of whole time evolution. 

After integrating the equation in t one obtains 


^—i{u!t+9) _ —2i{uit+9) _j_ ^ 

6 2 


Note 


x{t) 


xq + go au+{t) + 


(41) 


By equating —ih[u{t)x*(t) — u*{t)x{t)] computed by (39) and (41) to q{t) in (|T7|) , one 
obtains the following two equations: 

i«r-i/5r = 1, (42) 
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where xq denotes the constant shifting of annihilation operator, which gives a constant 
coherence, and {3 denotes the squeezing with respect to the true vacuum. 

The invariant annihilation operator is related to the annihilation and creation oper¬ 
ators at the true vacuum by the following relation: 


A 


i u*p — u* q -I- X* 
a*a + I3a\ 


(43) 




is the annihilation operator at g = 


where a = i u*^_{t)p{t) — u*^_{t){q{t) — go) + 
go with 0{x^) constant coherence. We cannot determine all the parameters without 
solving the full time evolution of the system, which needs numerical work. However, 
many important parts of the dynamics are determined without explicitly solving the 
full dynamics, e.g., the squeezing and coherence are determined just by the effective 
potential at t —>■ cxd, which is new results. 


4 Summary and discussions 

We have developed a LvN approach which gives the time-dependent Gaussian approx¬ 
imation. We, then, applied the method to a system with time-dependent slow rolling 
type potential, and found the time evolution of a Gaussian wave-packet initially centered 
at g = 0. We divided the evolution into three stages, and described solutions for each 
stages. We have calculated the period of the slow rolling phase which is signihcantly 
larger compared to the variationally determined value = \J^/ (Aa). We also displayed 
that many information of the hnal stage can be obtained without explicit solution of the 
full time evolution. Explicitly, we have shown that the invariant annihilation operator 
can be expressed as a linear combination of the creation and annihilation operators of 
true vacuum. We have also shown that the squeezing and constant coherence are de¬ 
termined only by the dynamics around the absolute minimum of the effective potential, 
independent of the details of the whole dynamics. 

This approach is the synthesis of the merits of the Gaussian approximation and LvN 
method. In the point of view of Gaussian approximation, it shows exactly the same 
potential in g and f as that of Ref. []^ , and in the point of view of LvN method, it was 
shown that the method can be applicable for general type of potential as a well posed 
approximation. There is, however, restrictions for the application range of the type of the 
potential. The Gaussian approach is not appropriate to describe the symmetry breaking 
potential from the start to the end, even though it is fairly accurate until f ~ This 

limitation shown in Ref. 0 was not cured in this paper. To overcome this limitation we 
need double Gaussian trial wave function or other approximations. 
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